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ABSTRACT 

We calculate corrections to the power spectrum predictions of weak lensing by large scale structure 
due to higher order effects in the gravitational potential. Using a perturbative approach to third order 
in transverse displacements, we calculate a second order correction to the angular power spectra of E 
and B mode shear and convergence resulting from dropping the so-called Born approximation, where 
one integrates along the unperturbed photon path. We also consider a correction to the power spectra 
from the coupling between lenses at different redshifts. Both effects generate B-modc shear and the 
latter also causes a net rotation of the background galaxy images. We show all these corrections are at 
least two orders of magnitude below the convergence or B-mode power and hence relevant only to future 
ultra high precision measurements. These analytical calculations are consistent with previous numerical 
estimates and validate the use of current large scale structure weak lensing predictions for cosmological 
studies and future use of B-modes as a monitor of systematic effects. 
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1. INTRODUCTION 

As experiments that measure the distortion induced in 
distant galaxy images by weak gravitational le nsing from 
the l arge-s cale structure of the universe (e.g., Blandforcj 
ct al| 1991; |Miralda-Escude| 1991; [Kaiseij 1992) move from 
the c urrent discovery phas e flBacon etaif 2000; [Kaiser et aj 
2000; |Wittman ct~a| 2000; [Van Waerbeke et al| 2000) into 
the precision measurement phase, it will become increas- 
ingly important to quantify and separate subtle cosmolog- 
ical, astrophysical and instrumental effects that alter the 
statistics of the lensing observables. 

Current predictions of weak lensing statistics and their 
utility in measuring fundamental cosmological parameters 
are based on several assumptions. In this Letter, we dis- 
cuss two of these assumptions and the extent of their va- 
lidity for calculations of the power spectra of the shear 
and convergence. In the so-called Born approximation, 
one integrates the lensing distortion over the unperturbed 
photon paths. We relax this assumption using the trans- 
verse deflection as a perturbative parameter. We show 
that the correction to the Born approximation results in a 
generation of curl, or magnetic-like, B-modes in the weak 
lensing shear. Absent in first order lensing contributions, 
B-modes are often used as a monitor of subtle systematic 
errors in the data. Indeed, data from the current gen- 
eration of surveys r outinely show B-modes in the shear 
field (e.g. Pen et al 2002). It will be important for future 
surveys to know the level at which lensing itself produces 
B-modes. 

An additional assumption is the neglect of the coupling 
between lenses at two different redshifts. Unlike the cor- 
rection to the Born approximation, the lens-lens coupling 
results in both B-mode generation and a net rotation of 
galaxy images. Since the rotational effect has been mea- 



sured in numerical simulations (Jain et a] 2000), it provides 
a check on the validity of our perturbative technique. 
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The physical mechanism behind the two corrections 
discussed in the present pap er were first con sidered by 
Bernardeau et al (1997) and Schneider et al (1998). In 



both these studies, these two second order corrections were 
discussed as a possible contribution to the three point 
statistics of weak lensing. The correction to the weak lens- 
ing convergence skewness, however, is below a few percent. 
In this paper, we concentrate on the two point statistics 
which requires a third order expansion in the potential to 
complete a second order expansion in the power. 

For illustrative purposes, we calculate these higher or- 
der effects for a flat ACDM cosmology throughout, with 
parameters Q c = 0.3, fib = 0.05, fiA = 0.65, h = 0.65, 
n = 1, Sh = 4.2 x 10 -5 and no tensor contribution. 

2. STATISTICAL PROPERTIES 

In this section, we define the statistical properties of the 
lensing observables in angular Fourier space. In general, 
we can write the weak lensing deformation matrix that 
maps between the source (S) and image (I) planes, 5xf — 
Sx j , as 



1 — k — 71 —72 — u> 
-72 + OJ 1 - K + 71 



(1) 



where all components are functions of position on the sky 
n. Here, k is the convergence and 7, are the two shear 
components and u>, the antisym metric component induces 
a rotatio n in the images (e.g., [Bartclmann & Schneider 



2000 and Mclliei 1999 for recent reviews), 
tion component may be isolated as 

VVb(n) = Sab - A ab (n) . 



The deforma- 



(2) 



Following Kaiser (1998), we define the Fourier representa- 
tion of the deformation under the flat-sky approximation 
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Statistical homogeneity requires that its two-point corre- 
lation satisfy 



(v>: 6 (i)<Mi')> = (2^) 2 £(i - v)c abcd (i) . 



(4) 



The various components of C a bcd define the two point 
statistical properties of the lensing observables. It is useful 
first to reexpress the shear components of ip a bO) m terms 
of t he coordin ate-f ree E and B mode r eprese ntation (see, 
e.g. |Stebbins| 1996; [Kamionkowski et al| 1998; |Hu fc White] 
2001) 
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cos2<^7i(l) 
cos 2^72(1) 



sin 2^72(1) . 
sin20;7i(l) . 
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where cos (f>i = I • Xi where Xi is one of the orthogonal 
directions on the sky. 

The power and cross spectra of the observable fields a, 
j3 = k, e, (3, u> are defined as 



<a*(W)> = (2tt) 2 <5(1 
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Since the power spectra are functions of the magnitude 
of l only, one can take l || Xi to simplify calculations of 
the power spectrum without loss of generality. With this 
restriction, — and we obtain for the power spectra 
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and similarly for the cross spectra. Cross spectra between 
k or 6 with (3 or lj vanish if the statistical properties are 
invariant under inversion of the coordinates. 

3. CALCULATIONAL TECHNIQUE 

It remains to evaluate the deformation power spectrum 
Cabcd to second order in perturbation theory. The de- 
formation tensor for sources at a sin gle redshift is given 
implicitly by (see e.g. Schneider et al 1998) 



i/>a&(n, Xs ) = 2 d X g(x, Xs)$ ac(x; x)[$cb + tpcb(n, X )] , 



(7) 

where commas represents spatial derivatives. Repeated 
indices are summed over the 2 transverse directions. The 
presence of ip in the integral reflects a foreground lens af- 
fecting the deformation from a more distant lens - or "lens- 
lens coupling" . The lensing efficiency is 



g(x',x) 



d A (x')d A {x-x') 

d A (x) 

0, x'>x- 



X < X; 



(8) 



Here the gravitational potential $ is evaluated at a de- 
flected position 2 

x(n, x) = nd A (x) + #x(n, x) , 

2 Formally d A (x)n —> d A (x)nj_ + x n \\ m an open universe but 
components parallel to the fiducial line of sight drop out in the Lim- 
ber approximation. 



Sx a (n, X ) = -2 J dx'gix'^^-^A^X 1 ), 



(9) 



where the deflections are confined to the transverse plane 
and x is the conformal distance or lookback time, from the 
observer, given by 



X{z) 



dz' 



Jo H(z>)> 

and the analogous angular diameter distance 



(10) 



d A ( X ) = H^n K 1/2 S mh(H n^x) 



1/2, 



(11) 

with Clfc = 1 — 51 t ot as the space curvature parameter. We 
will occasionally suppress the radial/temporal argument x 
where no confusion will arise. 

The familiar form of the lensing observables comes about 
by keeping only first order terms in the potential. In other 
words, they are calculated under the so-called Born ap- 
proximation where the potential is evaluated on the un- 
deflected path and with the neglect of lens-lens coupling. 
Higher order corrections can then be calculated by itera- 
tive correction of x and i\) a b in Eq. (@)- 

We shall see that in general the deformation may be 
expressed as a line-of-sight projection of a source field 



V>ab(n) = 2 / dxg(x,Xs)Sab(nd A ;x.) , 



(12) 



whose power sp ectrum may be evaluated under the Limber 
approximation (Kaiser 1992) from the power spectrum of 
the source field 

(S* ab (l; X )S cd (V; x')) = (2^) 2 <5(1 - l')8( X - x') 



xd A ( x r 6 Pabcd(l;x), (13) 



to be 



C abcd (\) = 4 / d X 



g 2 (x,Xs 



For the first order term S 
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where Pj,$(fc) is the spatial power spectrum of the poten- 
tial fluctuations evaluated under the no n-linear scaling rc- 
lations for the density power spectrum ( Peacock & Dodds 
1996). 

As discussed in the previous section, one may evaluate 
the power spectra of the observables under the constraint 
li = I and I2 = 0, from which one immediately obtains 
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and C," 11 



while C, 



0Pi 



0. Here the 



subscripts nn denote nth order term expansion in the po- 
tentials. The lack of P-modes in the signal, to first order 
in the perturbations, has been considered as a possible test 
of instrumental and astrophysical systematic effects. 

In general if S a b is symmetric under interchange of a and 
b, power spectra involving the rotation lo vanish. If P abc d 
is symmetric under the exchange of all indices then C? K = 



Ct 



cr. 

-in? 



Finally if P 1122 = P1221 then Cp* + C?» = 



Cf + C'i p . These relationships are useful for checking 
consistency in the higher order terms. 
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3.1. Born approximation 

Corrections to the Born approximation can be calcu- 
lated by Taylor expanding the potential in Eq. (R) to sec- 
ond order in the deflection of Eq. (|)J) 



<Z>(nd A + Sx) = <£>(nd A ) + 5x a <f>, a (nd A ) 
+ ^5x a 5x b $. ab (nd A ) - 



(17) 



The last term is effectively third order in the potential but 
must be kept since it can couple with the first order de- 
formation in the power spectrum. The corrections may be 
represented as second order and third order source fields 

^ 2 b ) (n) = -2$, afcc (nd A )B c (n), 
S$ (n) = 2<t>, abcd (nd A )B c (n)B d (n) , (18) 
with the Born correction 

B c (n- X )^ |d X , 5(x',x)^^$, c (nrf A (x , );x')- (19) 

The second-second order terms in the power spectrum 
become 



(22) = / fl_ 
abed I (- 27r )2 



I'Xl'Ml'-l'YMQ'J"), (20) 



with 1" = 1 — 1'. The mode coupling integrand is given by 



,9 2 (X,X') 



In terms of the lensing observables, 



X 



Cf 22 =4 / d X 



9 2 (x,Xs) f d 2 V ;/4c , (22)c (22) 



(21) 



d A { X f J (2nr ~° P 
x(l' -l") 2 M(l',l"). (22) 

The geometric factors for the specific observables are 
G( 22 )=cos^, 



M 22 ) ■ j. 
= sin (pv, 

G (22) = q 



(23) 



where we have used the condition </>; = 0. Notice that the 
parity violating cross terms vanish after integration over 
the azimuthal angle. 

Likewise the first-third order source term becomes 



P 



(13) 



abed 



~4:l a l b l c ld 



2;/ 



(1 • l') 2 M(l, I 1 ) . (24) 



Note that there is no mode coupling in the third order 
term. For the lensing observables, 



c\ 



■'hi . | / , u , 9 2 (x,Xs) f d 2 V ? 4 G (13) g (13) 



4 / dx- 



d A { X f J (2tt) 



x(M') 2 M(M') 



(25) 




Fig. 1. — Angular power spectra of weak lensing. The solid line is 
the well known first order result. The dot, dashed, long-dashed and 
dot-dashed lines are the total, Born approximation plus lens-lens 
coupling, second order contribution to lensing observables involving 
convergence, shear e-modes, shear /3-modes, and rotation, respec- 
tively. For comparison, the gray dashed line shows the contribution 
to shear e-modes from lens-lens coupling only; kk and /3/3 scale simi- 
larly. Note that the rotational power spectrum, lulu, is also generated 
by lens-lens coupling only. 



where 



= G (13) 



G 
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Notice that contributions to the power from modes with 
1' = I — 1" k, 1 in Eq. (22) is cancelled by contributions 
from Eq. (13) for kk and ee. Finite (3(3 power is generated 
by the Born corrections but is also highly suppressed geo- 
metrically. There is also no contribution to the rotational 
power from the second order Born approximation. The 
rotational power spectrum as measured in simulations by 
Jain et al (2000) should not be interpreted as a test of the 
Born approximation. 

3.2. Lens-Lens coupling 

Lens-lens coupling involves the iterative correction of 
the lensing deformation of distant lenses due to the defor- 
mation from foreground lenses in Eq. (|7|). The second and 
third order corrections become 

S%>(n) = -$, ac (nd A -,x) 

x2 / <frW>x)*,c6(ndA(x');X , )> 



S®(n) = $, ae (nd A ; X ) 

x2 / dx'g(x',x)<i>Mnd A (x');x') 



x2 J dx"g(x",x')*MnMx");x")- (27) 

The second-second order terms in the power spectrum 
become 



p(22) _ 4 
r abcd ~ ^ 



d 2 V 
(M 



i' a i' b %m'-\") 2 M{i',i") 



(28) 



with 1" =1—1' and recall that the mode coupling integrand 
was defined in Eq. (|2l]) . The first-third order term vanishes 
under the Limber approximation since lenses at the same 
redshift cannot lens each other. 

The results for the power spectra are 
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x(V -l") 2 M(l',l"), 
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with the geometric factors 

G< 2a > =C os(^ 
G( 22 )=cos(# + ^»), 



'/3 
/7( 22 ) 



■ Sin 0;/ 
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(30) 



Recall we have set <pi = 0. As derived, lens-lens coupling- 
generates power in all the lensing observables. Again, no- 
tice that the parity violating terms oscillate leaving negli- 
gible contribution after azimuthal integration. 



3.3. Born-Lens coupling 

There are also cross terms between corrections to the 
Born approximation and lens-lens coupling. The second- 
second order term follows directly from Eqs. ( O ) and (27) 



p(22) _ 4 



2// 



dH 



7 (i'Xix+mw-iyM(i',i") 



( 31 ) 

where again 1" = 1 — 1'. 

The second-second order terms for the lensing observ- 
ables are likewise simply related to the pure Born and 
lens-lens corrections 



Cf K22 = 4 



dx 



9 2 {x,Xs) 



d A (x) e 
x(l'-l") 2 M {I', I"), 



d* 1 ' 7/3;// r (22) r (22) 



(32) 



4. RESULTS & DISCUSSION 

Here, we presented a general discussion of the accuracy 
of first order calculations of statistical properties of weak 
lensing by large-scale structure. Using a perturbative ap- 
proach to third order in transverse displacements, we cal- 
culated the leading order corrections to angular power 
spectra of the shear and convergence that results from 
dropping the Born approximation. We also considered 
the coupling between lenses at two different redshifts in 
its contribution to the power spectra of lensing observ- 
ables. Both effects generate power in the B-mode shear 
which places an ultimate limit on its use as a monitor of 
systematics or for the search for weak lensing by gravita- 
tional waves. In addition, lens-lens coupling results in a 
net rotation of the background galaxy images. 

In Fig. |l|, we summarize our results and show the to- 
tal contribution to the lensing observables at the second 
order level in power from corrections involving the Born 
approximation and lens-lens couplings. As shown, these 
corrections are at least two orders of magnitude below the 
power in the convergence, or _B-mode shear, at the first 
order level. While these corrections are larger than cos- 
mic variance at I p> 100, they are unlikely to affect the 
interpretation of the next generation of surveys. Our an- 
alytical calculations are consistent with previous numer- 
ical estimates: the rotational power spectrum generated 
by the coupling between tw o lenses agree with numerical 
measurements by Jain et al (2000). These results also val- 
idate the use of current large scale structure weak lensing 
predictions for cosmological studies. The relatively small 
contribution to the lensing observables suggest that higher 
order corrections related to the Born approximation and 
coupling between two lenses are unlikely to affect the cur- 
rent estimates of weak lensing statistics and place only 
very mild limitations on the use of i?-modes as a monitor 
of systematics in future surveys. 
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where the geometric factors are simply products of the 
individual Born (B) and lens-lens (LL) factors 



^(22)^(22) _ r (22) 
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given in Eqns. (|23|) and (|30|). 
For the third order term 

S a f(n) = 2<I>, acd (nd A )B d (n) 



x2 / dx'g(x\x)$MndA(x'>,x'), (34) 



and 



p(13) _ 16 / fl_ 



(UblJ d + lJ'M(l-l') 2 M(l,l'). (35) 



This implies that all power spectra of lensing observables 
vanish after integration over the azimuthal angle </>c . 
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